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Supersymmetry, like Poincare symmetry, is softly broken at finite lattice spacing pro- 
vided the gaugino mass term is strongly suppressed. Domain wall fermions provide the 
mechanism for suppressing this term by approximately imposing chiral symmetry. We 
present the first numerical simulations of Af=l supersymmetric SU(2) Yang-Mills on the 
lattice in d=A dimensions using domain wall fermions. 

Supersymmetric (SUSY) field theories may play an important role in describing 
the physics beyond the Standard Model. Non-perturbative numerical studies of 
these theories could provide confirmation of existing analytical calculations and 
new insights on aspects of the theories not currently accessible to analytic methods. 
One such SUSY field theory that can be formulated and studied numerically on the 
lattice is the four dimensional Af—1 super Yang-Mills (SYM), which is just QCD 
with one flavor of adjoint Majorana fermions, called gluinos. In the traditional 
approach, the Wilson fermions are used to simulate the gluinos 0. Since Wilson 
fermions break chiral symmetry the gluinos will be massive, breaking SUSY, unless 
fine-tuned mass counterterms are introduced. Pioneering^ work using these methods 
has already produced very interesting numerical results 0'0. 

In this work, we summarize recent results of the first simulations using the lattice 
domain wall fermion formulation to represent the massless gluinos of SYM Bi. The 
idea behind domain wall fermions is to start with a formulation of massless fermions 
in one higher dimension and introduce domain walls so that chiral surface modes 
become exponentially bound to the walls. Pulling apart domain walls of opposite 
chirality, by increasing the size of the extra dimension, exponentially suppresses 
chiral symmetry breaking at a cost proportional to the size of the extra dimension. 
For reviews on DWF please see the LATTICE '00 review talk of Vranas and 
references therein. The possible use of DWF in SUSY theories has been discussed 
in earlier works Lffl and the methods used here are along these lines. Foil lists of 
references not included here for lack of space, please see the cited articles 

The Dirac operator in the adjoint representation of SU(iV) has an index 2Nv, 
where v is the winding of the gauge field. Classical instantons have integer winding 
and cause condensations of operators with 2N gluinos which anomalously breaks 
the U(l) -R-symmetry to %2N- For (xx) to condense, the remaining 7*2N symmetry 
must further break either spontaneously or anomalously to Z2. If the breaking is 
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anomalous, then the responsible gauge configurations must have fractional winding 
□. The existence of such gauge configurations has already been established It is 
our goal to distinguish between these two scenarios. 

All numerical simulations were performed on 8 4 and 4 4 Euclidean spacetime 
volumes with periodic boundary conditions using the inexact hybrid molecular dy- 
namics (HMD) R algorithm. The algorithm numerically integrates the classical 
equations of motion as part of generating a statistical ensemble with weights pro- 
portional to the fourth root of a two adjoint flavor Dirac determinant. For DWF, 
this weight is proportional to the weight for a single adjoint Majorana flavor. The 
integration step sizes were chosen such that systematic uncertainties due to numer- 
ical integration errors are negligible compared to statistical uncertainties. The 8 4 
volume simulations were run with /3=4/gr 2 =2.3, chosen as large as possible with- 
out entering the finite volume transition region. The 4 4 volume simulations were 
run with (3=2.1. Scaling arguments appropriate for weak coupling suggest that the 
lattice spacing at (3—2.1 is twice as large as at /3=2.3. 

To extrapolate the measured values of (xx) to the chiral limit, L s — >oo and 
toj— >0, simulations were performed in 8 4 volumes at fixed /3=2.3 while the size 
of the extra dimension L s was set to 12, 16, 20 or 24 and the bare mass rrif was 
set to 0.02, 0.04, 0.06 or 0.08. If the formation of the gluino condensate is due to 
spontaneous symmetry breaking, the lattice volume limits how small the dynamical 
qluino mass can be seijvithout losing the condensate: 12m e ff (xx) V ^> 1 (the 12 
is just normalization) tLJ. As m c s > nif, this limit is satisfied for all 8 4 simulations 
with m f > 0.02. 
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Figure 1: Extrapolated (xx) to 
L s —>oo limit vs. rrif and linear fit to 
rrif— >0 limit. 
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Figure 2: Dynamical (xx) a t m/=0 
vs. L s on 8 4 and 4 4 lattices. Curves 
are exponential fits. 



To estimate the gluino condensate in the chiral limit, we first extrapolate at fixed 
nif to the L s —>oo limit using the fit function (xx) = c o + c i exp(— C2i s )- The values 
of the extrapolated gluino condensate (with propagated errors) appear as points in 
figure |l|. These extrapolated values are further extraploted to the m/— >0 limit using 
a linear function (xx) \l ^oo = ^o+^i™/- The best fit function appears as the line in 
figure||, with bo—0- 00455(21). It is also reassuring to note that reversing the order 
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of limits, i. e. first m/— >0 at fixed L s then L s — >oo, yields a statistically consistent 
answer. 

Another approach to estimating the gluino condensate in the chiral limit is to 
actually perform dynamical simulations with m/=0. Since finite L s will induce an 
exponentially small breaking of chiral symmetry, the effective gluino mass will not 
be zero. However, the gluino mass should be too small to support spontaneous 
symmetry breaking. Additional simulations were run with L s set to 12, 16, 20 or 
24. The data are shown in figure ^. The curve is the best fit to an exponential 
fitting function with the extrapolated value of the condensate as shown. 

Surprisingly, both methods for estimating the gluino condensate produce con- 
sistent results within the statistical errors. Note that this is inconsistent with the 
notion of spontaneous symmetry breaking. Operationally, this result reinforces our 
claim that systematic uncertainties are still relatively small despite limited statis- 
tical precision. Further, it gives us some confidence that our fit functions are valid 
over the region of interest. 

To further check for spontaneous symmetry breaking of the Z4 symmetry, we 
measured (xx) 011 smaller 4 lattices with mj=0 and even larger values for L s . The 
data are shown in figure ^ with the best exponential fit and the extrapolated value 
for the condensate. This provides even stronger evidence that spontaneous symme- 
try breaking is not responsible for the formation of a gluino condensate, at least 
finite volumes. On these lattices 12m/V(xx) < 1; so analytical considerations 
suggest the support of (xx) must come primarily from topological sectors with 
fractional winding of v = ±1/2. 

The spectrum of the theory is of great interest but it was not possible to measure 
it on the small lattices considered here. Also, the gluino condensate was measured 
at only two different lattice spacings so extrapolation to the continuum limit to 
compare with analytical results is not possible. Future work could explore these 
very interesting topics. 
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